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We study the finite-size scaling of heavy-light mesons in the static limit. The most relevant effects 
are due to the pseudo-Goldstone boson cloud. In the HMChPT framework we compute two-point 
functions of left current densitities as well as pseudoscalar densitites for the cases in which some 
or all of them lay in the epsilon-regime. As expected, finite volume dependence turns out to 
be significant in this regime and can be predicted in the effective theory in terms of the infinite- 
volume low-energy couplings. These results might be relevant for extraction of heavy-light meson 
properties from lattice simulations. 
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1. Introduction 

The simulations of heavy-light (hi) mesons made out of a heavy quark (charm or bottom) and 
a light one (up, down or strange) on the lattice are challenging because they require very large 
volumes and small lattice spacings in order to keep systematic errors under control. 
Requiring that a < l/(2m^) and M//L > 4, where Mu is the mass of a pseudo-Goldstone boson 
made of two light quarks, would imply a simulation at a < 0.02 fm and L > 6 fm if Mu ~ M n , with 
an L/a > 300. 

If the heavy quark mass is sufficiently large a good effective description is provided by heavy quark 
effective theory (HQET) [1-3], which is a systematic expansion in the limit of infinite heavy quark 
mass. In this limit, the scale M/,/ (the mass of a heavy-light meson) disappears from the problem 
and the UV cutoff can in principle be as low as the cutoff used to describe light meson dynamics. 
Indeed this approximation has been extensively used to simulate heavy-light mesons in lattice QCD 
(see [4] for a review on heavy flavour phenomenology on Lattice QCD). 

On the other side, both if the heavy quark is treated in the static limit or not, an obvious question is 
whether we can do better concerning the constraint on the box-size. After all the finite-size scaling 
of heavy-light systems should be dominated by light pions physics, since these are the lightest 
modes in QCD. To the extent that pion physics can be described by chiral perturbation theory 
(ChPT), it is conceivable that finite-size scaling of heavy-light systems can be accurately predicted 
using ChPT, as the finite-size scaling of light mesons is [5, 6]. 

Having integrated the bottom out, we can take the lattice spacing as big as a < l/(2m c ) ~ 0.05 
fm; then in a lattice of length L = 2fm we can reasonably expect that ChPT still provides a good 
description of QCD. At this volume, a lattice computation with masses in the epsilon regime would 
require L/a > 40. 

In [7] we have computed the finite-size scaling (FSS) of heavy-light systems when the lightest 
pions are light compared to the inverse box size in two limiting situations. In the first one the 
heavy quark is significantly heavier than the light one, but still treatable in ChPT: this correspond 
to considering hi mesons in the mixed-regime introduced in [8] ; in the second one the heavy quark 
is static and therefore chiral dynamics can be treated in Heavy Meson Chiral Perturbation Theory 
(HMChPT): this correspond to considering hi static mesons in the £-regime. 
Even though these two situations are physically very different, the pion dynamics responsible for 
the finite-size scaling properties should be pretty much the same. It is therefore interesting to see 
explicitly how a quantitative matching of the fmite-size effects takes place, by comparing the finite 
volume dependence of correlation functions in ChPT and HMChPT. 

We present the result obtained for the left-handed current correlator to next-to-leading order in both 
effective theories and discuss its applicability to extract LECs from lattice simulations. Full results 
and details are presented in [7]. 

2. Heavy-light mesons in the mixed-regime of ChPT 

We consider N/, heavy quarks of mass nif, and Nj light quarks of mass m; with nij <C but 
such that both can still be treated in the context of ChPT. In practice we apply the power counting: 

m/~£ 4 , m h r^e 2 , L 1 ~^~£, (2.1) 
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that has been named mixed-regime in [8], to the common ChPT Lagrangian. While we refer to 
[8] for further details on the implementation of ChPT in the mixed-regime, both in the full and 
partially-quenched theories, we just remind that in this regime the zero modes of the light pions 
have to be treated non perturbatively. Following the same method we have computed the two-point 
correlation functions of two heavy-light left-handed currents to relative tff(e 2 ) order: 

Tr[T a T b ]Cj(t) = Jd 3 x(j$(x)J b (0)) (2.2) 

where J" = \j/T a y^P^\j/ and T a is a traceless generator of SU (Ni +Nh) with one index in the light 
and one in the heavy subsector. 



2.1 Results in the mixed regime of ChPT 

The result is: 



C J (t) = F ^Mf A) P(t,Mf A) )-^. 
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where we have defined the renormalized and M?n as: 



Nh ~ ^) k 00 (M 2 ,2M 2 ,t) + (Jj- + 1) k 00 (M 2 ,M 2 ,t) I, (2.3) 



F 2 ~ \ fa - ^) G (0,21$ - f G(0,0) - N -^G (0M) 
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ft, \6A// 4 4 

and we have defined the shorthands ft = m ; £V , M| = an d M 2 = N + ' N M\. Here we have used: 
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G(x,M 2 ) = -Y 
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N t 



EE{X) V^ (p 2 ) 2 F(p) 2N 2 VM 2 > F{P) p 2 + 2M 2 + p 2 
The time dependence is expressed in terms of the function: 

cosh[M(f -\t\)] 



(2.6) 
(2.7) 



P(t,M 2 ) = / d 3 xG(x,M 2 ) 



and its derivatives: 



2Msinh [^f] 

k Q0 (M 2 ,M 2 ,t) ee l -^i 2 ^(t,M 2 lp )^M P )- (p(tM P )^(t,M 2 2p ) + (Mi «-M 2 )) 1, 



(2.8) 



(2-9) 
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where M ap = \/M 2 + p 2 . In this regime we also need: 

(Pit M 2 )M 2 \ 
MM 2 ,M 2 V)+ 2T m2 1 j ' ( 2 - 10 ) 

The result in Eq. (2.3) represents the FSS of kaon-like states (m/, = m s and m/ = m u = m^) for 
various situations: 

• 2 + 1 dynamical simulations: Nh = l,Ni = 2; 

• partially quenched (PQ) simulations where the h quarks are quenched and the / quarks are 
dynamical by taking the replica limit Nf, — > of Eq. (2.3); 

• PQ simulations where the / quarks are all quenched or PQ, while the h quarks are dynamical. 
In this case, the appropriate value of 2V/ must be taken, but also the zero-mode integrals 
((U Q + U%) u) need to be properly defined [9, 10] 1 . 



3. Heavy-light mesons in the HMChPT 

The effects of pion dynamics in the properties of static heavy-light mesons can be predicted in 
HMChPT [12-14]. In [15], chiral corrections to B parameters in the ^-regime were studied. As 
far as we know, the e -regime has not been explored yet. However the technology we have used was 
developed in [16] to perform £ -regime calculations in baryon ChPT 

The most relevant complication is that in the HMChPT Lagrangian half integer powers of the 
pseudo-Goldstone field appear, so that the observables are not just functional derivatives of the 
zero modes partition functional (which can be calculated explicitely). We dealt with this problem 
by working with an explicitely parametrization for the Goldstone bosons manifold. We restricted 
to the Ni = 2 case. The result can be written in the form 2 : 



a 2 6(t) 



N,=2 



exp 



-am( £ », 



1 + 



3 1 

4 JFL)'- 



(H(t,L,T)+glH'(t,L,T)) 



(3.1) 



where 



H(t,L,T) = L 2 ( L hl (L\ + * £ P ( f ,p2) -G(0,0) 



1 



H'(t,L,T) = - £(p(,, p 2)-P(0,p 2 )) 



\t\ 1 

T ~ 2 



1 

12 



(3.2) 



and 



3g 2 , 
8F 2 L 3 ' 



a = Fp^j2M P = V2F P */y/Mp^ 



(3.3) 



'Note that one cannot consider a fully quenched theory with N/, = Ni = on the basis of Eq. (2.3), because the 
singlet has been integrated out[l 1]. 

2 Note that we use Vjit) for HMChPT and Cjif) for ChPT. The reason for this will be clear in section 5. 
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The coupling g n is responsible for the scattering processes BB*K and B*B*K. It has been computed 
on the lattice by several authors; the most recent determination has been done in [17] where the 
value g K = 0.44 has been found. It is interesting to stress that the dependence on m/ and on the 
topological sector disappears after all diagrams are summed. 

This result may be used to analyze the scaling for correlators of B or D systems. To analyze the 
B s and the B* , or to take into account the effects due to the strange quark explicitely for B or B* 
systems, one has to add a quark in the p-regime, as we discuss in [7]. However in these systems 
the influence of the scalar and axial resonances is not negligible, as argued in [18], and the validity 
of HMChPT without including them is questionable. 

In Fig. 1 (left), we show the ratio of the finite- volume to infinite volume correlator in Eq. (3.1) at 
t = 1 fm as a function of L for two different geometries and two values of g K . Corrections are 
0(4%) at 2 fm, and the depedence on g K is mild. In Fig. 1 (right), we show the time dependence 
of Eq. (3.1) after factoring out the exp(— AM^ e h) (we will see later that in any real fit, AM would 
renormalize the hi meson mass). 



mi ~ e , N, = 2, N s = 



m, ~ e 4 , N, = 2, N s = 
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Figure 1: Left: ratio of ffj(t = Ifm) at fixed volume normalized to the °° volume result as a function of L 
for two lattices with T = L (solid) and T = 2L (dashed), and for g K ~ 0.44 (thick lines)[17] and g K = (thin 
lines) and typical values of E = (250MeV) 3 ,F = 90MeV. Right: 4%(t)/(a 2 exp(-AM £ t)) as a function of 
t/T for T = L = 2 fm, and for g K = 0.44 (thick line) and g K = (thin line). 



4. Matching 

Dominant finite-size effects in QCD are due to pion dynamics, since these are the lightest 
degrees of freedom. It is therefore expected that the finite-size scaling of heavy-light systems does 
not depend on the large energy scales related to the heavy quarks, ie. Mi, or . This must be the 
case as long as those scales are significantly larger than L _1 . Whether these scales are much larger 
also than the QCD scale so that the static limit (HQET) is a good approximation, or not, should not 
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matter a priori for the finite-size scaling properties, because the volume dependence arises from the 
propagation of the light degrees of freedom. 

The leading finite volume effects are therefore expected to come from the fact that the heavy meson 
can emit and absorb a pion. The probability for this to happen can however depend on the heavy 
mass scale. Close enough to the chiral limit, the masses of pseudoscalar mesons are suppressed by 
the spontaneous breaking of chiral symmetry, so, for example, we do not need to include the vector 
mesons in the effective theory, because they are much heavier and decouple. To the contrary, in 
the limit m/j —> °° pseudoscalar and vector mesons are degenerate because the interaction between 
quark and antiquark inside the meson becomes spin independent, so they both need to be considered 
in HMChPT. The presence of heavy-light vector resonances can modify the finite volume effects 
indirectly by inducing unsuppressed contributions to pion/heavy-light meson scattering. We will 
see that indeed the finite-size corrections in HMChPT and mixed ChPT match up to corrections 
proportional to g\. 

We show now how the matching works. Given any meson two-point function, the first point to 
realize is that a finite static limit is recovered after factorizing out the leading e~ Mt , where M is the 
mass of the hi meson. What can be matched is the dependence of the correlators on the volume, that 
is L and T and the masses of the up, down and strange quarks because these are explicit degrees of 
freedom in both effective theories. Moreover, since we only consider the static limit of HMChPT, 
we have to drop from the ChPT result those contributions that are suppressed by negative powers 
of nth. In order to match eqs. (3.1) and (2.3), the L, T dependence must be the same in both cases. 
For the mixed ChPT framework, we split the contribution due to the heavy pions from the rest in 
Eqs. (2.4) and (2.5) and write: 

N, 



if A) =F K,iV ( )-- iVr- 6(0,0) + %'), (4.1) 



Mf A) = M 2 h (m h , Ni) + 0(m h l ) . (4.2) 

F and Mf, have absorbed the dependence on the heavy quark mass. These expressions just remind 
us that F and £ depend on the number of quarks that are considered in the effective theory. The 
static limit of the mixed ChPT case in Eq. (2.3) is, for t > 0: 



Cj(t) F 2 M h 



1 + 2W(^) 



(4.3) 



exp(-M h t) 4 

where H(t,L, T) is the function of Eq. (3.2). 

For Nj = 2, the result is identical to the NLO prediction % (t) in HMChPT (Eq. (3.1)) with the 
following identifications: 



a = F P ^2M P = F y 2M/,, g n = 0. (4.4) 

The fact that at NLO we have to put g n = to match the two expressions reflects the fact that the 
vector meson is integrated out in the chiral theory. 

5. Finite-size scaling of heavy-light mesons in lattice QCD 

As we have seen above the matching of finite-size effects of heavy-light correlators in HM- 
ChPT and ChPT works as expected. We would be interested however to use these results to predict 
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the finite-size scaling of these correlators computed in lattice QCD. We can do this including a 
relativistic heavy quark or in the static limit. In both cases we expect that for sufficiently large time 
separation: 

cT{t)=^ a ^V^))lat^%{t) x J_exp(-M ? ), (5.1) 

where M is the lightest heavy-light meson mass in the case of a relativistic heavy quark or the 
so-called static energy, E stat = M/,/ — m/, in the lattice static limit. 

Note that the value of E stat is not predicted by HMChPT, however in general we can write: 

E stat = E^ t +AM, (5.2) 

where E^} t is the value the static energy would have in the chiral limit, while AM contains the 
chiral corrections that are predicted by HMChPT. In the £-regime case they are given in Eq. (3.3). 
In practice this means that to fit a correlator using Eq. (3.1) one has to determine four parameters: 
a, F, E^lf and g K . It remains to be seen what the stability of such fits is in practice. 
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